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The alternating and non-alternating harmonic sums and other algebraic objects of the same equivalence class are 
connected by algebraic relations which are induced by the product of these quantities and which depend on their 
index class rather than on their value. We show how to find a basis of the associated algebra. The length of the 
basis I is found to be < 1 /d, where d is the depth of the sums considered and is given by the 2nd Witt formula. 
It can be also determined counting the Lyndon words of the respective index set. There are two further classes 
of relations: structural relations between NiELSEN-type integrals and relations due to the specific structure of 
Feynman diagrams which lead to a considerable reduction of the set of basic functions. The relations derived can 
be used to simplify results of higher order calculations in QED and QCD. We also report on results calculating 
the 16th non-singlet moment of unpolarized structure functions at 3-loop order in the MS scheme. 



1. INTRODUCTION 

Multiple nested alternating and non-alternating 
harmonic sums Sai,...,a„iN) [1~3] emerge in pcr- 
turbative higher order calculations within QED 
and QCD for massless fermions, 
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Here, are positive or negative integers and 
A'^ is a positive even or odd integer depend- 
ing on the observable under consideration. One 
calls n the depth and J2t=i Wk\ the weight of 
a harmonic sum. Harmonic sums are associ- 
ated to Mellin transforms of real functions or 
SCHWARTZ-distributions f{x) e S'[0, 1] [4] 

Sa,,...,aAN)= f dx fa,....,aA^) (2) 

Jo 

which emerge in field theoretic calculations. Fi- 
nite harmonic sums are related to harmonic 
polylogarithms Hh^,,,,^h^{x) [5]. Their 1/(1 ± 



a:)-weighted Mellin transform yields harmonic 
sums. The inverse Mellin transform relates 
the harmonic sums to functions of Nielsen in- 
tegrals [6] of the variable x at least for all sums 
of weight w < 4 as shown in [3], and associated 
generalizations for higher weight. Nielsen inte- 
grals are a generalization of the usual polylog- 
arithms [7]. In the limit iV — > c» the conver- 
gent multiple harmonic sums, i.e. those where 
ai 7^ 1, yield (multiple) Zeta-values Cai,...,a„, 
which are also called Euler-Zagier sums [8] . A 
generalization of both harmonic polylogarithms 
and the Euler-Zagier sums are the nested Z- 
sums [9], which form a HoPF algebra [10,11] and 
are related to GONCHAROv's multiple polyloga- 
rithms [12]. For a recent review see [13]. 

Higher order calculations in massless field the- 
ories are either performed in Mellin-A^ space 
referring to harmonic sums or in the space of the 
momentum fractions x representing the results in 
terms of NiELSEN-type integrals. The principal 
complexity is determined by the amount of pos- 
sible terms contributing. In the case of the 2- 
loop coefficient functions in momentum fraction 
space [14] 77 different functions occurred, cf. [3]. 
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This number compares to the amount of all pos- 
sible different nested harmonic sums up to weight 
w = 4, 80 = S'" — 1. For the 3-loop anomalous 
dimensions [15] one expects the contribution of 
a wide class of the w = 5 harmonic sums and 
for the 3-loop coefficient functions of the w = 6 
harmonic sums, which means 242 or 728 sums, re- 
spectively. These sums are not independent but 
are connected by different kind of relations. In the 
present paper we summarize a first class of rela- 
tions recently being discussed in Ref. [16], the so- 
called algebraic relations. It turns out that these 
relations emerge from the index-structure and the 
multiplication relation of the objects considered 
and are widely independent of other properties of 
the harmonic sums. In this way an equivalence 
class of even more objects is defined having the 
same properties or can be found as special cases 
thereof. One example is the set of the harmonic 
polylogarithms [5]. 

To obtain manageable expressions it is of im- 
portance to apply all these relations through 
which the number of basic functions to be re- 
ferred to is considerably reduced. A further re- 
duction of the number of basic functions follows 
due to structural relations [17] and the symmetry 
of Feynman amplitudes. 

Experience shows that the Mellin space rep- 
resentation yields simpler expressions in gen- 
eral [18], which would not be seen easily working 
in X space. To obtain as simple as possible ex- 
pressions is of special importance because of the 
fact that data-analyses require compact results 
in a:-space, which can be obtained using analytic 
continuations for the basic sums [19] and per- 
forming a single numeric Mellin inversion [20] 
for the whole problem. Since the evaluation of 
precise analytic continuations needs special effort 
any possible reduction carried out before is of 
help. 

2. RELATIONS BETWEEN HARMO- 
NIC SUMS 

There are three types of relations between har- 
monic sums. 

1. Algebraic relations due to the shuffle alge- 
bra of harmonic sums under their multipli- 



cation, cf. [16,21,22]. 

2. Relations implied by structural relations of 
NiELSEN-type integrals [17], which are re- 
lated to harmonic sums via a Mellin trans- 
forms after being weighted with the denom- 
inators 1/(1 ± x) 

3. Symmetry relations encoded in Feyn- 
man diagrams for the respective quanti- 
ties, which cancel a series of harmonic sums 
present in individual diagrams via algebraic 
relations. 

The product of two finite harmonic sums (1) 
yields 



^ sign(ai)'i ^ 
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-Sai,...,an (h) Sb2,...,bm (^2) 



[sign(ai)sign(6i)] 



■■•.an(0'S'62,...,6„(0 ■ 

(3) 

We introduce the shuffle product LLI of a single 
and a general finite harmonic sum 

+Sbi,ai,b2,---,bm{^) + • • • + 'S'6i,62,...,6m,ai(-^) 

(4) 

which is a linear combination of the sums of depth 
TO + 1 . The shuffle product of two harmonic sums 
of depth n and m, Sai,...,a„{^) and Sb^^,,,^b^{N), 
is the sum of all harmonic sums of depth m + n 
in the index set of which occurs left of aj for 
i < j and likewise for bk and bi for /c < I. As 
an example the shuflile product of two threefold 
harmonic sums is given by 

,0.4, ,a^,aQ 

(N) 

+Sa (TV) 
{N)+Sa (TV) 

~\~Sai ,a4 ,a2 ,0.5 ,ae ,^3 {-^) ~^ ,04 .05 ,05 .02 .as (-^) 

•Jai ,a4,a5 ,02,0-3 ,ae (TV) 
{N) + 'S'o4,a5,ai,a6,a2,a3 

(TV) 

H~'S'o4,a5 ,ai ,a2 ,a6 ,03 

(TV) 



3 



(N) 

~\~Sa4,ai ,a5,a2,a3,a6 

{N) 

,ai,a2,a5,03,a6 

(iV) 
(5) 

Finally one establishes a system of linear equa- 
tions in which the linear elements of the shuffle 
produets form the variables and a polynomial out 
of harmonic sums of lower depth forms the in- 
homogeneity. We furthermore consider all index 
permutations. This system contains all algebraic 
relations. In Ref. [16] all solutions for harmonic 
sums up to depth d = 6 were given. This complies 
to the level of sophistication needed to reduce the 
corresponding relations which emerge for mass- 
less 3-loop coefficient functions. 

3. NUMBER OF ALGEBRAICALLY IN- 
DEPENDENT HARMONIC SUMS 

Let us consider the index set of a harmonic sum 
of depth d. One may consider this set word 
w or a, non-commutative product of letters of an 
ordered alphabet 21 = {a,b,c,d, . . .}. Any word 
can be decomposed into three parts 



w = pxs 



(6) 



a prefix p, a suffix s, and the remainder part 
X. Among all words w the Lyndon words, cf. 
e.g. [23], are those being smaller than any of its 
suffixes. 

According to a Theorem by Radford [24] the 
shuffle algebra discussed above is freely generated 
by the Lyndon words, i.e. the length of its basis 
is given by the number of Lyndon words. We 
would like to count the number of Lyndon words 
for index sets, where the same letters can emerge 
repeatedly. The corresponding relation is due to 
Witt [25] and will be called 2nd Witt formula 



ln{ni, ...yUq) 
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with n = J2k=i^k- Here iJ,{d) denotes the 
MoBius function. One may derive ln{ni, ■ ■ ■ ,lq) 



using the generating functional 
1 

. . — -r„ 

n=l 



1 — Xi 



oo / ^ 



i„(n.) 



Z^fc=l 



dh 



Note that (7) is related to the Gauss-Witt re- 
lation mentioned by Hoffman [21] for the num- 
ber of basic multiple Zeta-values of weight w for 
Vr7.,; > if all cases for fixed weight arc summed 
over. An even more strict relation in the inclu- 
sive case has been conjectured by Zagier [8] and 
Broadhurst and Kreimer [26] in the case of 
multiple Zeta-values and verified up to w = 12. 

Let us come back to Eq. (7). We can draw some 
immediate conclusions out of this relation. If the 
numbers tii have no common divisor larger than 1, 
the number of the basis elements compared to the 
number of all objects equals l/rf, where d denotes 
the depth of the index set. In case of common 
divisors larger than 1 we checked that the basis 
is always shorter for all depths up to d = 10, 
sec [16]. 



Weight 


# Sums 


# Basic Sums 


Fraction 


1 


2 





0.0 


2 


8 


1 


0.1250 


3 


26 


7 


0.2692 


4 


80 


23 


0.2875 


5 


242 


69 


0.2851 


6 


728 


183 


0.2513 



4. OTHER RELATIONS 

Further to the algebraic relations structural re- 
lations between the generating functions of har- 
monic sums, the 1/(1 ± .x) weighted Nielsen type 
integrals, exist. These relations reduce the num- 
ber of basic harmonic sums further [17]. As an 
example we mention the relation 



1U|(^ 
2 1 - a;2 



Li2(a;) Li2(-a;) 



+ 



1-x 
Li2(x) 



+ 



+ 



1 - 

Li2(- 



X 

-x) 



l + x 



(8) 



1-1- a; 

The structure of Fcynman diagrams is further- 
more selective on the type of harmonic sums 
which may occur for the final physical results. 
As mentioned before, a;-space calculations of the 
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two-loop Wilson coefficients resulted in the oc- 
currence of nearly as many different functions as 
combinatorially possible. Studying the Mellin 
transforms of these functions in detail [18,27] one 
finds that a series of harmonic sums occurs with a 
partial index symmetry such that certain classes 
of basic functions do not contribute. In the case 
of two loop Wilson coefficients the function 



ln(l 



1 + x 



(9) 



which is mathematically irreducible, does not 
contribute, although it is of weight 2 and func- 
tions up to weight 4 span the basis. Likewise 
other irreducible lower-weight functions, which 
cannot be eliminated by structural relations do 
not occur. The function (9) is, however, instru- 
mental to build the basis spanning the set of func- 
tions required for the 3-loop anomalous dimen- 
sions. 

5. A QUADRATIC LAW ? 

The final analysis of the anomalous dimensions 
and Wilson coefficients [18, 27, 28] shows that 
one may come to the following representation : 

1. Single harmonic sums S±k{N) and their 
analytic continuation oc ■ijj''^^{z) are re- 
garded as well-known functions. Likewise 
any derivative of the analytic continuation 
of the Mellin transforms over [0, 1] of the 
basic functions introduced subsequently are 
considered to be trivial. 

2. The analytic continuation of the Mellin 

transforms of the basic functions are mero- 
morphic functions with poles at the inte- 
gers bound from above by a positive num- 
ber. They can be represented by factorial 

series [29]. 

3. The one-loop anomalous dimensions re- 
quire w=l harmonic sums, likewise the 1- 
loop Wilson coefficient are expressed by 
at most w=2 harmonic sums, which are all 
trivial in the above sense. 

4. The two-loop anomalous dimension are rep- 
resented as polynomials of trivial functions 



but the Mellin transform of one non- 
trivial function of up to w=3 



Li2(a;) 
1 + x'' 



[!]• 



(10) 



5. The two-loop Wilson coefficients are rep- 
resented as polynomials of the Mellin 
transform of the functions mentioned and 
the four additional basic functions 



Li2(a:) Li3(.-r) Si,2(x) 

1 — X ' 1 + X ' l±x ^ 



[18, 27]. (11) 



6. The three-loop anomalous dimensions are 
represented as polynomials of the Mellin 
transform of the functions mentioned and 
the nine additional basic functions 

ln(l -|- x) Li4(a;) 



1 + x 
Si,3(a:) 
1 + x 



l±x' 

1±X ' 



S2,2{-x)-Ul{-x)/2 Ul{a 



l±x 



1 + x ' 



[28]. (12) 



In this way the number of basic functions con- 
tributing to the single scale quantities discussed 
up to weight w = 5 leads to the following reduc- 
tion from the combinatorial complexity 



2-3 



e{w -2)[w- 2Y 



(13) 



cf. [28], i.e. the number of basic elements 
grows only quadratically, rather than exponen- 
tially. This relation will be examined for the case 
of three-loop coefficient functions, i.e. for w=6, 
very soon. 

6. THE 16th MOMENT OF THE NON- 
SINGLET THREE LOOP ANOMA- 
LOUS DIMENSION OF Fi{x,Q'^) 

Previous fixed-moment calculations of anoma- 
lous dimensions and WiLSON coefficients ex- 
tended to the moments = 14 for the non- 
singlet structure function ^1,1,(0;, Q^) and A'' = 
13 for xF3{x,Q^) [15]. Weeks before the cal- 
culation of the complete result for the 3-loop 
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non-singlet anomalous dimensions [31] was fin- 
ished ^ the calculation of the 16th moment of the 
non singlet anomalous dimension of the structure 
function Fi(x,Q'^) and the WiLSON coefficients 
for F2^l{x, Q^) was started using the MINCER algo- 
rithm [30]. The result of this calculation provides 
an independent test of the complete results. The 
calculation was performed using several power- 
PC's partially linked to a 4.2 Tbyte raid system. ^ 
The calculational time for both the projections 
g^v and P^^Py of the forward Compton ampli- 
tude amounted to 564 CPU days. The 16th mo- 
ment of the non-singlet anomalous dimension for 
PP(a;,Q2) reads, cf. [32]: 



+ 



64419601 
6126120 
117652: 



Cf 



(14) 



•:5810303 



112588038763200 
21546159166129889 



CfNf 
Cf Ca 



484994628518400 
3689024452928781382877 
" 459818557352009856000 



Cf^(15) 



(2) _ / 59290512768143 
'•'i® ~ I 1563722760600 

58552930270652300886778705063429867 \ 
3451337970612452534317096673280000 j 

■Cf' 

( 15018421824060388659436559 
^ 579371382263532418560000 

/ 1670423728083984207878825467 
\^ 6488959481351563087872000 

59290512768143 ^ \ ^ 2 
+ 3127445521200 ^-'j^^^^ 

5559466349834573157251 2 
2069183508084044352000 " " 



corresponding publication on the 3-loop WiLSON coef- 
ficients will appear later. 

^Wc thank S. Wicsand and P. Wcgncr for technical assis- 
tance and U. Gensch and C. Spiering for their support of 
this calculation. 



1229794646000775781127856064477 
30335885575318557435801600000 

59290512768143 \ 2^ 
1042481840400 ^^j^'' 

( 71543599677985155342551355451 



+ 



\^ 938967886855098206346240000 

64419601 



765765 



Cs Cf'Nf 



(16) 



and agrees with the result in [31] for A'' = 16. 
The corresponding moments of the WiLSON coef- 
ficients C2,l{x,Q^) are given in [32]. 

7. CONCLUSIONS 

The product of finite alternating or non 
alternating harmonic sums is given by the shuf- 
fle product of harmonic sums and polynomials of 
harmonic sums of lower depth. These representa- 
tions imply algebraic relations between the har- 
monic sums. If one considers all harmonic sums 
associated to an index set {ai, . . . , Ofc} one may 
express these sums by a number of basic sums. It 
turns out that this number is given by the 2nd 
Witt formula which counts the number of Lyn- 
don words corresponding to the respective index 
set. The set of these Lyndon words generates in 
this sense all harmonic sums of this class freely. 
By solving the corresponding linear equations we 
derived the explicit representation of all harmonic 
sums up to depth d — 6 without specifying the 
indices numerically and gave all expression which 
are structurally needed to express the sums up to 
weight w = 6. The counting relations for the ba- 
sis of the finite harmonic sums were given up to 
depth d= 10. The relations derived hold likewise 
for other mathematical objects obeying the same 
multiplication relation or a simpler one, which is 
being contained, as that for harmonic polyloga- 
rithms. This is clue to the fact that the relations 
derived depend on the index set and the multi- 
plication relation but on no further properties of 
the objects considered. 

The ratio of the number of basic sums for a 
given index set and the number of all sums is 
mainly determined by the depth d rather than 
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the weight of the respective sums, due to the pre- 
factor in the WiTT formula. Modifications 
occur due to common non-trivial divisors of the 
numbers of individual indices in the set being con- 
sidered. Up to (i = 10 we showed that the frac- 
tion of basic sums is always <l/d compared to all 
sums. The use of these algebraic relations leads 
to a considerable reduction in the set of functions 
needed to express the results of higher order cal- 
culations in massless QED and QCD and related 
subjects. Further reductions are due to structural 
relations and symmetries in the set of Feynman 
amplitudes calculated. For practical applications 
such as the description of the QCD scaling vio- 
lation of the structure functions in deeply inelas- 
tic scattering the harmonic sums occurring in the 
Mellin N space calculation have to be translated 
to a;-space by the inverse Mellin transform. For 
this reason the respective harmonic sums have 
to be analytically continued in the argument N 
to complex values, which requires a high effort 
using numerical procedures. It is therefore re- 
commended to use as many as possible relations 
between the N space objects before the last step 
is being performed. 

The 16th moment for the parity-conserving 
non singlet structure function was calculated. It 
provides a test on the complete result, which has 
been derived recently for the anomalous dimen- 
sions, and on the upcoming result for the coeffi- 
cient functions. 

Acknowledgment. This paper was supported 
in part by DFG Sonderforschungsbereich Tran- 
sregio 9, Computergestiitztc Thcorctische Physik 
and EU grant HPRN-CT-2000-00149. I would 
like to thank S. Moch, V. Ravindran, T. Riemann 
and J. Vermaseren for discussions. 



REFERENCES 

1. A. Gonzalez- Arroyo. C. Lopez, and F.J. Yn- 
durain, Nucl. Phys. B153 (1979) 161; 

A. Gonzalez-Arroyo and C. Lopez, Nucl. 
Phys. B166 (1980) 429. 

2. J. A.M. Vermaseren, Int. J. Mod. Phys. A14 
(1999) 2037. 

3. J. Bliimlein and S. Kurth, DESY 97-160, 



hep-ph/9708388; Phys. Rev. D60 (1999) 
014018. 

4. L. Schwartz, Theorie des Distributions (Her- 
mann and Cie, Paris, 1951), Vols. I and II; 
V.S. Vladimirov, Gleichungen der Mathematis- 
chen Physik, (DVW, Berlin, 1972); 

K. Yoshida, Functional Analysis (Springer, 
Berlin, 1978). 

5. E. Remiddi and J. Vermaseren, Int. J. Mod. 
Phys. A15 (2000) 725. 

6. N. Nielsen, Nova Acta Leopold. XC, 121 
(1909); 

S. Kolbig, Siam J. Math. Anal. 17 (1986) 
1232. 

7. L. Lewin, Dilogarithms and Associated Func- 
tions (Mac-Donald, London, 1958); Polyloga- 
rithms and Associated Functions, (North Hol- 
land, New York, 1981); 

A. Devoto and D.W. Duke, Riv. Nuovo Cim. 
7 N6 (1984) 1. 

8. L. Eulcr, Novi Comm. Acad. Sci. Petropol. 
20 (1775) 140; 

D. Zagier, First European Congress of Math- 
ematics, (Birkhauser, Boston, 1994), Vol.11, 
p. 497. 

9. S. Moch, P. Uwer, and S. Weinzierl, J. Math. 
Phys. 43 (2002) 3363; 

S. Weinzierl, Comput. Phys. Commun. 145 
(2002) 357. 

10. H. Hopf, Ann. of Math. 42 (1941) 22; 

J. Milner and J. Moore, Ann. of Math. 81 
(1965) 211; 

M.E. Swcedler, Hopf Algebras, (Benjamin, 
New York, 1969). 

11. C. Kassel, Quantum Groups, (Springer, Berlin, 
1995). 

12. A.B. Goncharov, Math. Res. Lett. 5 (1998) 
497, 

http : / / www . math . uiuc . edu/K-Theory/ 0297. 

13. M. Waldschmidt, J. Theor. Nombres Bor- 
deaux 12 (2000) 581; Multiple polyloga- 
rithms, Lectures at the Institute of Math. Sci- 
ences (Chennai, Nov. 2000); updated version 
http : //www . math . jussieu . f r/~miw/ 
articles/ps/mpl .ps; 

V.V. Zudilin, Uspekhi Mat. Nauk (Russian 
Math. Surveys) 58 (2003) 3. 

14. W.L van Neerven and E.B. Zijlstra, Nucl. 



Phys. B382 (1992) 11. 

15. S.A. Larin, T. van Ritbergen, and J. A.M. 
Vermascrcn, Nucl. Phys. B427 (1994) 41; 
S.A. Larin, P. Nogueira, T. van Ritbergen, 
and J. A.M. Vermaseren Nucl. Phys. B492 

(1997) 338; 

A. Retey and J. A.M. Vermaseren, Nucl. Phys. 
B604 (2001) 281. 

16. J. Bliimlcin, Comput. Phys. Commun. 159 
(2004) 19. 

17. J. Bliimlein, DESY 04-064. 

18. J. Bliimlein and S. Moch, in preparation. 

19. J. Bliimlein, Comput. Phys. Commun. 133 
(2000) 76. 

20. J. Bliimlcin and A. Vogt, Phys. Rev. D57 

(1998) 1; D58 (1998) 014020. 

21. M.E. Hoffman, J. Algebra, 194 (1997) 477. 

22. M.E. Hoffman, math. QA/0406589 and these 
Proceedings; Pacific J. Math. 152 (1992) 
275; J. Algebraic Combin. 11 (2000) 49; 
math.qA/9907173; math. NT/0401319; 
Ngoc Minh Hoang, these Proceedings and ref- 
erences therein. 

23. C. Rcutcnaucr, Free Lie Algebras, (Calendron 
Press, Oxford, 1993). 

24. D.E. Radford, J. Algebra, 58 (1979) 432. 

25. E. Witt, Journ. Reine und Angcw. Mathc- 
matik, 177 (1937) 152; Math. Zeitschr. 64 
(1956) 195. 

26. D.J. Broadhurst and D. Kreimer, Phys. Lett. 
B393 (1997) 403. 

27. J. Bliimlein and V. Ravindran, these Proceed- 
ings and in propartion. 

28. J. Bliimlein and S. Moch in preparation. 

29. N. Nielsen, Handbuch der Theorie der Gamma- 
funktion, (Tcubncr, Leipzig, 1906); 

E. Landau, Uber die Grundlagen der Theorie 
der Fakultatenreihen, S.-Ber. math.-maturw. 
Kl. Bayerische Akad. Wiss. Miinchen, 36 
(1906) 151. 

30. S. A. Larin, F. V. Tkachev, and J. A. M. Ver- 
maseren, NIKHEF-H-91-18. 

31. S. Moch, J. Vermaseren and A. Vogt, Nucl. 
Phys. B688 (2004) 101. 

32. J. Bliimlein and J. Vermaseren, in prepara- 
tion. 



